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A ROBUST INVERSION METHOD FOR QUANTITATIVE 3D 
SHAPE RECONSTRUCTION FROM COAXIAL EDDY-CURRENT 

MEASUREMENTS 

HOUSSEM HADDAR*, ZIXIAN JIANG*, AND MOHAMED KAMEL RIAHI* 


Abstract. This work is motivated by the monitoring of conductive clogging deposits in steam 
generator at the level of support plates. One would like to use monoaxial coils measurements to 
obtain estimates on the clogging volume. We propose a 3D shape optimization technique based on 
simplified parametrization of the geometry adapted to the measurement nature and resolution. The 
direct problem is modeled by the eddy current approximation of time-harmonic Maxwell’s equations 
in the low frequency regime. A potential formulation is adopted in order to easily handle the complex 
topology of the industrial problem setting. We first characterize the shape derivatives of the deposit 
impedance signal using an adjoint field technique. For the inversion procedure, the direct and 
adjoint problems have to be solved for each coil vertical position which is excessively time and 
memory consuming. To overcome this difficulty, we propose and discuss a steepest descent method 
based on a fixed and invariant triangulation. Numerical experiments are presented to illustrate the 
convergence and the efficiency of the method. 

Key words. Electromagnetism, non-destructive testing, time harmonic eddy-current, Inverse 
problem, shape optimization. 
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1. Introduction. Non-destructive testing using eddy-current low frequency ex¬ 
citation are widely practiced to detect magnetite deposits in steam generators (SG) 
in nuclear power plants. These deposits, due to magnetite particles contained in the 
cooling water, usually accumulate around the quatrefoil support plates (SP) and thus 
clog the water traffic lane. Many methods and softwares based on signal processing has 
been developed in order to detect deposits using standard bobbin coils and are widely 
operational in the nuclear industries (see for instance the database of nondestructive 
testing m and references therein). Estimates of the bulk amount of deposits enable 
to supplement a chemical cleaning process, which in some cases may be ineffective 
where it leaves significant deposits in the bottom area of the SP foils. The presence of 
such deposits generates a reduction and re-distribution of the water in SG circulation 
and can cause flow-induced vibration instability risks. This may harm the safety of 
the nuclear power plant. 



Fig. 1. Three dimensional mesh of the SG and the SP clogging : coils (pink), tube (blue), SP 
(grey) and deposit (red). 
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3D shape reconstruction from eddy-current measurements 


In order to obtain better characterizations than those provided by model free 
methods, we present and discuss a robust inversion algorithm (for non destructive 
evaluation using eddy current signals) based on shape optimization techniques and 
adapted parametrizations for the deposit shapes. An overview of techniques for non 
destructive evaluations using eddy currents can be found in [7] and we also refer 
to [23] . [20 and [26, for further engineering considerations. For other model based 
inversion methods related to eddy-currents we may refer, without being exhaustive, 
to 0 El m H3 E]- In the medical context, several inverse source problems related 
to eddy-current models have been addressed: non-invasive applications for electroen¬ 
cephalography, magnetoencephalography 0 (see also pQ) and magnetic induction to¬ 
mography hmd. 

Stated more precisely, the inverse shape problem we shall investigate aims at 
retrieving the support of a conductive deposits using monostatic measurements of 
coaxial coils and their computable shape derivatives. Our work can be seen as an 
extension of m to a realistic 3D industrial configuration. Although the deposit 
geometry can be an arbitrary three dimensional domain, the available (monostatic) 
measurements can only give qualitative information on the width. Since the objective 
is to detect the possibility of clogging at the support plate, we found it appropriate 
to consider a deposit concentrated in only one of the opening regions at the support 
plate (See Figs. [l]). The geometrical parameters are then the deposit width at (at 
most) one measurement position. In practice, it turned out that a relative robustness 
with respect to noise can be achieved if one shape parameter correspond with two 
vertical positions of the coils. In order to speed up the inversion procedure, we are led 
to consider a fixed geometrical mesh (adapted to the chosen parametrization). This 
allows us to obtain an inversion procedure which is not very sensitive to the number of 
measurements. Moreover, in order to avoid troubles due to changes in the conductive 
region topology, we adopted a vector potential formulation of the 3D eddy-current 
model. A careful study of the shape derivative of the solution to this formulation is 
conducted. For related shape derivatives associated with Maxwell’s equations we refer 
to [9j [15] . We here treat the potential formulation of the eddy current problem. This 
derivative allows us to rigorously define the adjoint state, needed to cheaply compute 
the coils impedances shape derivatives. 

The geometrical setting of the industrial configuration is depicted in Fig. [1] We 
denote by Q the computational domain, which will be a sufficiently large simply 
connected cylinder. It contains a conductor domain flc composed of the tube, the 
support plate and eventually a deposit on the exterior part of the tube: U 

where t stands for the tube, d for the deposit and p for the SP. The insulator 
domain Q \ Vtc is split into two parts: that indicates the region inside the tube 

where the coil (thus the source J) is located and Q v that denotes the insulator outer 
region (where the deposit can be formed). For our purpose we introduce the surface 
T = dfld H dVt v that denotes the interface between the deposit and the insulator. 

Let us now briefly describe the 3D eddy-current model, which derives from the 
full Maxwell’s equation in the time harmonic low frequency case and the adopted 
formulation of this problem. Given the bounded domain Q C M 3 , we recall the time- 
harmonic Maxwell equations: 

curl H + (icje — cr)E = J in D, 

curlE — icj/iH =0 in D, 

and on the boundary dfl we impose a magnetic boundary condition Hxn = 0, where 
n stands for the outward normal to the boundary dd. Here H and E denotes the 
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magnetic and electric fields, respectively. J is the applied current density, e is the 
electric permittivity, /i is the magnetic permeability and cr is the electric conductivity. 
In our case, the applied current density has support strictly included in the insulator 
fit (interior of the tube). By neglecting the displacement current term, we formally 
obtain the eddy-current model, which reads: 

curlH = crE + J in £2, 
uH =—curlE in £2, 

1 UJ 

We refer to the monograph [2j for an extensive overview of eddy-current models and 
formulations. In this paper we adopt a potential formulation in which we look for the 
magnetic vector potential A and electric scalar potential V (only defined on Qq) that 
satisfies 


( 1 . 2 ) 


/iH 

curl H = crE 

E 

div A =0 


n x 




curl A 

in f}, 

CrE + J 

in fi, 

iu:A + W 

in Q c , 

0 

in Q. 

0 

on 

0 

on 9S2, 


where the equation ( 1 . 2)4 stands for the Coulomb gauge condition. The boundary 
condition JH5 stands for the magnetic boundary condition and the boundary con¬ 
dition (1.2)6 is equivalent to eE • n = 0. The electric scalar potential V is determined 
up to an additive constant in each connected-component of VLq , which has a connected 
boundary. 

Notice that from Maxwell-Ampere equation ( |l.l| )i we get 


(1.3) 


curl (fi 1 curl A) — cr(iu;A + VF) = J in Q. 


In the following, the space il(curl;S2) indicates the set of real or complex valued 
functions v G (L 2 (S2)) 3 such that curlv G (L 2 (0)) 3 and define 

:= {v G il(curl, fi), div v = 0 in f2, v • n = 0 on 

For a vector magnetic potential A G A(Q), an electric scalar potential V G H 1 (Qc) / 
C (the quotient by constants is relative to each connected component separately) and 
a test function G X(Q) the weak formulation of (1.3) reads 


(1.4) 


/ “ 

Jn V 


— curl A • curl T' dx 
n V 


- f a (iwA • ^ + W • *) dx = f J • 'i>. 

Jn c Jsi 


Moreover, for any test function <3? G H 1 (Qc) /C the weak formulation of the necessary 
condition —div (crE) = divJ writes — f Qc crE • V^ds = f Qc J • V^d s. Therefore 
using (1.2 ) 3 we obtain: 

(1.5) 


[ (io;A + W) • V0dx = [ 
J Oc J O 


= / J • dx. 


Following |2 Chp-6] (and references therein), by introducing a constant /i*, repre¬ 
senting a suitable average of fi in Q, the Coulomb gauge condition ( 1 . 2)3 can be 
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incorporated in equation (1.4) in the following way 

( 1 . 6 ) 

/ — curl A-curl T' dx-\ -/ div Adiv dx— / a (icoA • T' + W • T') dx = / J\l> 

Jn V Jo Jn c Jo 


and the variational space A'(fi) would then be replaced by := If (curl, fi) D 

ifo(div,fi) or equivalently by if 1 (Si) 3 since the domain Si is convex and sufficiently 
regular. Indeed, div A = 0 is verified in the weak sense. Combining equations ( |1.6| ) 
with we can obtain a symmetric variational formulation as follows 

(1.7) S(A,V;¥,$) = / J-^ds- — [ J-V$d;r V($,J)eQ, 

Jo la; JOc 

where the sesquilinear form <S is defined by: 


S(A, V; Tq 4>) := f (— curl A • curl T' + —div Adiv \Eu dx 

( 1 . 8 ) ' n i >J f _ 

+ — / cr(io;A + W) • (icj\l/ + V4>) dx. 

1UJ Jo c 

The coercivity of S on H 1 ^) 3 x if 1 (Sic) /C (see for instance [2j Chp-6]) ensures the 
well-posedness of the problem. 

This paper is organized as follows: after this introduction, we state in Section [2] 
the nonlinear shape optimization problem by the introduction of the misfit function, 
which depends on the shape of the defect and in particular its eddy-current signal 
response. We derive, in Section [3j the adjoint problem which is based on the shape 
derivative of the misfit function. At the end of this Section we explicitly formulate 
the shape gradient via the adjoint problem. In Section [4j we present and explain the 
algorithm of steepest descent based on the use of fixed predefined grid. With Section[5j 
we conclude the paper with numerical experiments that illustrate the robustness of 
the method. Some technical materials related to shape derivative are reported in the 
appendix for the readers’ convenience. 

We conclude this section with the introduction of some useful notations. We 
denote by [•] the jump across the interface T: [F] = lim^o F{x + tn) — lim^o F(x — 
tn ) \/x £ T, we recall here that n denote the normal to T pointing outside Qd- 

For any vector A and differentiable scalar V, we respectively denote the tangential 
component and the tangential gradient on some boundary or interface having a normal 
n by A r := A — (A • n)n and \7 r V := W — d n V ■ n. We finally shall use the notation 
G(ft) := U(Sl) x H 1 (Q c ) /C. 

2. Statement of the inverse problem. 

2.1. Impedance measurements. The deposit probing is an operation of scan 
with two coils introduced inside the tube along its axis from a vertical position ( min to 
a vertical position £ max . At each position ( £ [Cmin, Cmax], we measure the impedance 
signal Z((). According to [7, (10a)], in the full Maxwell’s system, the impedance 
measured in the coil k when the electromagnetic field is induced by the coil l writes 

A Zki = YY \2 fdn d (^° x H k — Efc x H^) • ndS ;where and are respectively 

the electric field and the magnetic field in the deposit-free case with corresponding 
permeability and conductivity distributions /i°, cr°, while E&, H& are those in the 
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case with deposits. Using the divergence theorem, we obtain the following volume 
representation of the impedances 

A Z kl » f div (E? x H fc - E fc x H?) dx 

= T7i2 f ( curl E ° ' H fe - E ? • curl H fc - curl E fc • H? + E fe • curl H?) dx 

\ J \ Jn d 

(2.1) “ ~Tm I (i~ - ^o)curlE fc -curlE° - iw(<r - <r°)E fc • E°) dx. 

1 CJ I J I Jn d V M M / 

In the last equality we used the eddy-current model CD- Furthermore, using the 
relation E = ic^A + VU we replace the electric field E by the vector potential A and 
we thus obtain the following shape dependent impedance measurement formula^ 

( 2 . 2 ) 

A Zkiifld) 

= [ ((- - -q) curl A/e • curl A° - ^(cr - cr°)( i(jA k + VT4) • (icjA° + VU/ 3 )^ dx. 

\ J \ Jn d V V M J 

2.2. A least squares formulation. Let us denote by £) the impedance 

response signals of a probed deposit that we would like to estimate. We shall 
use the shape dependent form in the impedance signal response in order to convert 
the signal anomaly to a shape perturbation. This inverse problem will be solved by 
minimizing a least square misfit function representing the error between computed 
and observed signals integrated over the coil positions. This misfit function is defined 
as follows: 


(2.3) 


r Cm; 

m d )= / 

J Cm i n 


12(0*0 - C)| 2 dC, 


where Z is either Zfa or Zp 3 according to the measurement mode used in practice: 
ZFA(^d) ■■= J(AZn(n d ) + AZ 21 (O d )), or Z F3 (n d ) := i(A Z u {Q d ) - AZ 22 (O d )). 
Minimizing this functional using a steepest descent method requires a characterization 
of its derivative with respect to perturbations of 0^. This is the objective of next 
section. 

3. Adjoint problem and explicit formulation of the shape gradient. We 

shall first study the shape derivative of the solution (A, V) with respect to deforma¬ 
tions of the deposit shape. This derivative will then allow us to obtain an expression of 
the cost-functional derivative. A computable version of this derivative is then derived 
through the introduction of an adjoint state. 

3.1. A preliminary result on the material derivative. In this part, we 
formally derive the expression of the material derivative of the solution to the eddy- 
current model on a regular open set with constant physical coefficients /i, a. This 
result will be used in next sections to obtain the material derivative of the eddy- 
current model with piecewise constant coefficients as well as the shape derivative 


* Let’s recall the fact that cr° is an e-conductivity. 
E° = icjA° + W°. 


Hence the electric field E® has a sense with 




6 


3D shape reconstruction from eddy-current measurements 


of the impedance measurements. We begin by introducing the shape and material 
derivatives m Section 6.3.3]. For any regular open set ft C M 3 , we consider a 
domain deformation as a perturbation of the identity Id + 6 : ft ftg,x y , 

where 0 G C := (C 2 (M 3 ; M 3 )) 3 is a small perturbation of the domain. To make a 
difference between the differential operators before and after the variable substitution, 
we denote by curl^, div^, V x the curl, divergence and gradient operators on ft with 
x-coordinates, and respectively by curl y , div y , V y those on fig with //-coordinates. 
For any (A (fie), V(Clo)) defined on fie, we set 

A C url 0 9 ) := (I + V0)*A(n fl ) o (Id + 0), 

A div (0) := det(/ + V6>)(/ + VO^A^g) o (Id + 0), 

V v (0) := v(n e )o(u + e). 


These quantities conserve the corresponding differential operators in the following 
sense (see for example [22J (3.75), Corollary 3.58, Lemma 3.59]) 


(3-1) 


I + VO 

det (I + VO) 

1 

det (I + VO) 

(I + VOr^VviO) 


curia, A cur }(0) 

div cc A div (0) 


= (curl y A(fl 0 )) o (Id + 0), 

= (div yA(fle)) o (Id + 0), 
= (V y y(^))o(Id + 0), 


where V0 := (— —)ij is the Jacobian matrix. 

OQCj 

In order to simplify the notation we use curl, div and V for respectively curl*,, div a, 
and Va;. 

Let (A (fl),V(fl)) be some shape-dependent functions that belong to some Banach 
space W(O), and 0 G C a shape perturbation. The material derivatives (B(0),t/(0)) 
of (A, T7), if they exist, are defined as 


(3 2) J ^" cui3 (0) -A-curl (0) + B(0) + o(0 ) — A(Q) + B(0) + o(0\ 

\u v (0) =Vv{0) + U(0) + o{0) = V{n) + U{0) + o{0), 

We also define the shape derivatives (A'(0), V'(0)) of (A, V ) by 

3 f A \0) := B(0) - (0 • V)A(fi) - (V0)*A(O), 

lU'(0) :=U(0)-O-VV(Q). 


The derivative Bdi v (0) of A which conserve the divergence operator is given by 
(3.4) B div (0) := B (0) + (div 01-VO- (V0)*) A(fi). 


Using the chain rule, in any open set of fl fl Q,g we formally have 

(3.5) A(fie) = A(f2) + A'(0) + o(0), 

(3.6) A div (0) = A(fi) + B div (0) + o(0), 

(3.7) V(fl e ) = V(Si) + V'(0) + o(0). 


To ease further discussions, in particular the derivation of the variational formulation 
(3.18) from (3.15), we give a preliminary result. Assume that the coefficients /u and 






H. Haddar, Zixian Jiang and M.-K RIAHI 


7 


a are constant on Q. We set a shape-dependent form 


(3 ' 8) 

A(n)(A,V;^,$) := / — curl A • curlSf dx -— f + VV) • (icj\l/ + V<1>) dx. 

Jn L Jn 


Compared to the variational form S defined in 
the penalization term ) -1 div Adiv dx. 

Lemma 3.1. Let Lt be a regular open set, fi > 0 and a > 0 constant on 
Ll and Id + 0 : Ll Lie a given deformation. Let (A,V) = (A(Q),V(Q)) and 
(\F,<F) = (\F(H), <f>(H)) be some shape-dependent functions with sufficient regularity. 
We assume that the material derivatives ( B(6 ), U(0)) of (A, V), the shape derivatives 
(A'(0),V f (0)) of (A, V) and the material derivatives (rf(0),xW) °f defined 

with ( |3.2| ) exist. If (A(£l),V{Q)) satisfy in the weak sense 


(1.7), the above form A(11) get rid of 


(3.9) 


curl(/i -1 curl A) — a(iuA + VV) = 0 in H, 

< div A = 0 inQ, 

^affujA + VV) ■ n = 0 on dll, 


then the shape derivative of A(Ll) that we denote by A'(H)(0), i.e. A(Qq)(A, V] \F, <F) = 

A(n)(A,y ; ^,^) +A / (n)(^)(A,y ; ^,^) +o(0), satisfies 

(3.10) 

A'(fi)(0)(A,K;¥,$) = A(fi)(A'(0),r(0);^,$) +A(fi)(A,F;77(0),x(0)) 

+ / — (0 • curlA)(n • curl\I>) ds 

Jan L 

— / cj{ti • 9){\ujA t H- V T y) • (icj\F r -f- V r <F) ds. 

1UJ Jan 


The proof of this Lemma is given in the Appendix. 


3.2. Material derivative of the solution to the eddy-current problem. In 

this part, we show the existence of the material derivative of the solution to the eddy- 
current problem with respect to a domain variation, and give its weak formulation with 
a right hand side in the form of some boundary integrals. We rewrite the variational 


formulation of the eddy-current model (1.7) on Qq. For any test functions (\F, $) G Q 

(3.H) 

s(A(n e ),v(n e y,^(n e ),me)) = 


/,( 


— curl y A • curl^ \l> -f-div y Adiv ) dy 


I 


\L 

-f 

^ Jn 


M* 


r(icjA ■ 


J • \l> dy - 

1 UJ 


VV) • (L;^ - 
J • V4> dy. 


V<f>) dy 


^ce 


We choose the test functions as follows (so that their material derivatives vanish) 
= (I -b VO) t 'S'(Qo) o (Id + 0).<L> = <£(£1(9) o (Id + 0). Since the supports of J and 6 
are disjoint, i.e. supp(J) D supp(0) = 0, the right-hand side of the weak formulation 
(3.11) writes simply: 


[ J • dx — ~~ [ 
Jn lu; Jn c 


(3.12) 


J • V<L> dx. 
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We consider the following term which conserves the divergence operator 


*div (0) := det(/ + V0)(/ + V0)“ 1 (/ + V0)“ t O' 

= det(/ + V0)(I + V0) _1 ¥(fie) o (Id + 0), 
div^div (0) = det(J + V0)(div v ®(n e )) o(Id + 0). 


By variable substitution y = (Id + 0)x, the left hand side of (3.11) can be written as 
(3.13) 

f /1 (I + V0)«(J + V0) ,^1_i_ 

7nU | det(7 + V0)| curl + ^ | det(/ + V0)| 

+ — / o-|det(/ + V0)|(/ + V0)- 1 (J + V0)- t (iwA curl +W v ) • (i^TW)da;. 
1W 7n c 


div A div div ^ 


div 


dx 


Theorem 3.2. Let 0 £ C a domain perturbation. Let fi > 0, a > 0 belong 
to L°°(Q). We recall that J E L 2 (Q) 3 has compact support in Ll s C Ll v and satisfies 
div J = 0 in Ll s and supp( J)Dsupp(^) =0. If (A(Q),V(Ll)) = (A cur i(0), Vv(0)) is the 
solution to the eddy-current problem ( |1.7| ) and (A(Qq),V(LIq)) = ((/+V0) _t A cur \(6)o 
(Id + 0) _1 , Vs/ { 6 ) o (Id + 0) _1 ) the solution to the problem flim} ; then 


lim 

ll^llc^O 


(Acurl(^) 


^-curl(O), Vs/{0) ~ Vy(0)) 


Q 


= 0 . 


Proof. We recall that A(fi) and A (Qq) satisfy the Coulomb gauge condition on 
Ll and on Qq respectively: div A(fi) = 0 on div^A (Lie) — 0 on Qq. From the weak 
formulations (1.7), (3.11) the identities ( |3.12| ), (3.13) and the developments in (A.l) 
we obtain 
(3.14) 

s(A curl (0) - A curl (0), V v (0) - W(0); IP,#) 

7 l 


+— 


(div 01 — V0 — (V0)*) curl A cur i(0) • curllZ' dx 


- [ ■ 
Jn c 


>'(—div 01 + V0 + (V0) t )(iwA cur i(0) + W v (0)) • (iwlP + V<P) dx + o(0). 


Obviously the right hand side of the above equality goes to zero as ||0||c —> 0. Since the 
form S is coercive (see [U Section 6.1.2]), this implies || (A cur i(0) — A cur i(0), Vy ( 0 ) — 
by(0))|| s —^ 0 as J]011c —> 0. □ 

Theorem 3.3. Under the same assumptions as in Theorem \ 3.2\ the material 
derivative of the solution (A(Q),V(£l)) to the eddy-current problem ( |1.7~| ) with respect 
to a domain variation Id + 0 exists. If it is denoted by (B(0),U(0)), then 


1 


lim 

e\\c^o \ 0 


(A curl (0) - A curl (0) - B(0), V v (0) - V v (0) - U(0)) 


= 0. 


Proof. Let (B(0), U(0)) the unique solution in Q to the weak formulation 


(3.15) 


<S(B(0), 17(0); 0,<P) = L(0,0) V(«P,d>) e Q, 


























H. Haddar, Zixian Jiang and M.-K RIAHI 


9 


where 

(3.16) 


L(¥, g>) := f n - (div 01 - V6» - (V0)‘) curl A • curl O' dx 
C / 

1 -div ((div 01 — VO — (V0)*)A) div&dx 


/ 

Jn 


n d 0* 


+ — / <r(-div0/ +V0 + (V0) t )(iwA + W) • (iwS'+ V<2>) dy. 

lw ^O c 

Let Bdiv (0) defined by (3.4). Then we can rewrite the weak formulation (3.15) as 


L 


— curlB(0) • curl#' + — divB div (0)div\P ) dx 
o C* 

1 


— J n cr(kjB(0) + W(0)) • (iu;^ + V<£)dx 


(3.17) 


1CJ 


/ 

JO 


1 


+ 


— (div 01 — V0 — (V0)*) curl A • curl \P dx 
~ o C 

J Qc cr(—div 0/ + V0 + (V0)*) (iw A + W) • (i^TW) dx. 


1CJ 


From (3.6) and the Coulomb gauge conditions satisfied by A(fi) and A(Q$) we deduce 
that divBdiv (0) = o(0). Considering the fact that (A, V) = (A cur i(0), Vy(0)), ( |3.14 ) 
and ( 3.17| ) yield 

s(A cml {0) - A curl (0) - B (O),V v (0) ~ V v (0) - U(0 ); <2>) + o(0) 

= [ l(div07 — V# — (V0)*)(curlA cur i(0) — curlA cur i(0)) • curl^drc 
JO M 

+ 4 f cr(—div 01 + V6» + (V0)‘) fiw(A curi (6>) - A curi (0)) + (VVv(6>) - W v (0)) • (k^TW) dx. 
Jo c V 7 

Theorem 3.2 implies that the right hand side of the above equality is of order o(0) as 
||0||c 0. The coercivity of S ensures the result as stated. □ 

Proposition 3.4. Under the same assumptions as in Theorem \3.2\ we assume 
in addition that fa, a are piecewise constant and constant in each subdomain (Q s , Ut t , 

Qd> ^v or Qp)- If the domain perturbation 0 has support only on a vicinity of the 
interface T between the deposit domain Qd and the vacuum fl v (T = D dUt v ) and 
vanishes in £} s , then the material derivatives (B(0),U(0)) of (A,V) satisfies 


S(B(0), £7(0); #C #>) = £(#C <£) V(#C <P) G Q, 


(3.18) 

where 

(3.19) 

■= 

[ (- curl((0 • V)A + (VO) 1 A) • curliF + —div ((0 ■ V)A + (V0) f A)div<F ) dx 
Jn d VC C* 


+ — [ cr( iw((0 • V)A + (VO)* A) + V(0 • W)) • (iw«P + V<£) da: 
1CJ Jo c V / 

+ / — (0 • n)(n • curl A)(n • curl#') ds 

Jr L/C 

4 /(O-n)M(i 

iCJ Jp 


+— 


(icjA r + V r P) • (icjlPy + V r <£) ds. 
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Proof. Let A := {s,£,d, v,p} a set of indices with its elements indicating the 
different sub-domains as well as the corresponding permeabilities and conductivities. 
We rewrite left-hand-side of the variational formulation (3.11) as 

S( A(fie), V(Q e y, A i(^)( A » V-, ¥, $) 

ieA 

+ [ —div y A(Sle) * div^^) dy. 

According to the definition of the test functions <h(U< 9 )), their respective 

material derivatives vanish. Since (A(Ufl), V(Qq)) satisfy both ( |1.6| ), we can apply 
Lemma 3.1 to the terms A^(O^), which yields the shape derivative 

]Ta'(^)(0)(a,V^) 

ieA 

= ]T A^OW), 17(0); 0) + Y, A«(a«)(-(0 • V)A - (V0)*A, -(0 • W);<F,$) 


ieA 


ieA 


L 


— {6 • curl A) (n • curllZ') 

/i 


ds — — f (O' n)[cr](i 
\id Jy 


\UjA t + V r V) • (ico^r + V r <2>) ds 


(3.20) 


= <S(B(0),f/(0);«F,<£) - f —divB(0)div?Fdx 

Jo 

+ £ Ai(ni)(-(0 • V)A - (V0)*A, -(0 • W); <F, «P) 


ieA 


L 


(0 ' n)(n • curl A) (n • curllZ') ds 
— f (0 • n) [a] (iwA T + V T V) • (iwO> T + V T <?) ds. 

iCJ Jp 


In the last equality we have used the transmission conditions [n • curl A] = [n x 
(/i _1 curlA x n)] =0 on T. Using the identities ( |A.1| ) and the Coulomb gauge 
condition div A = 0, one verifies that on each subdomain (i G A) of U 


(3.21) 


div ((div 0/ - V0 - (V0)*) A) = -div ((0 • V)A + (V0)‘A). 


From the derivation of L(>F. <I>) (3.16) and the equality ( |3.21[ |. one easily deduces that 

the shape derivative of the penalization term [ Q —div y A(Q^j • div y SI/ (fl#) dy is 

6 fA* 

[ — divB(0)div^da; + [ —div ((div07 - V0 - (V0)*)A Jdiv^dx 
Jn M* Jn d M* V / 

(3.22) = [ — divB(0)div^drc- [ 

Jn M* Jn c 


1 


—div ((0 ■ V)A + (V^) t A)divlZ >r dx. 


We easily get from (3.20) and (3.22) the variational formulation (3.18) with £(\P,<I>) 
given by (3.19). □ 
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3.3. Expression of the impedance shape derivative using the adjoint 
state. Now we shall give a new expression of the impedance measurements using 
the above results and the adjoint state. We recall the expression of the impedance 
measurements (2.2) 


A Z k i(£ld) = Tj^ 


f- 7 T) curl A*- • curl A? 


id a>' 

- 2 - (<X - (7°) (iwAfe + Wife) • (iwA? + W,°) j Ax. 

VjJ J 


Proposition 3.5. Let (Ak,Vk) be the solution to the variational formulation 
<H3 with coefficients (i, cr, and (A®, VJ°) the solution to with coefficients pf, a 0 
which do not depend on the deposit domain Lid- Let (A' k ,V k ) be the shape derivatives 
of (A, V ). Under the same assumptions as in Proposition ^. 4\ the shape derivative of 
the impedance measurement A Z k i(flf) is given by 

(3.23) 

A Z' kl {Q d )(0) 

= A f (f - To ) curl A k • curl A °i -fy- a°)& A 'k + vv fe) • ( iujA ° + VV; 0 )) Ax 
\ J \ Jvt d \ L L 1Ld J 


1C d 

w 


— curl A k • curlA° — ^[a\(iu;A kr + V r T4) • (iu;A^ r + V r Vj°) J ds. 

Jr V IPi 1UJ J 


Proof From (2.2) one has 

|J| 2 


-AZi 


ki 


(n d ) = A(n„) (A fc , v k -, a?, -v; 0 ) - A 0 (n d ) (a?, v?-,A k , -v k ), 


where A and A 0 are the forms defined in ( |3.8| ) with respectively the coefficients (/i, a) 
and (/i°,cr°). As (A kl 14) (resp. (A °, Vft )) satisfies (3.9) with constant coefficients 
(/i, a) (resp. (/i°, cr 0 )) in Lemma [3d] implies 

(3.24) ff-A Z' kl (n d )(9) 

1 UJ 

= A'(00(0) (Afc, 14; A?, -if) - /k' 0 (Q.d){0) (a?, Vf A~ k ,-Vk) 

= A(fi d ) (A ' fc (0), V£(0);A°, -Vf) + A(Q d ) (A fc , V*; B°(0), -tf(0)) 

- A 0 (n d ) (a°>), V) o, (0); Ai, -R) - A 0 (fi d ) (a°, VJ°; Bffl), -R(0)) 

+ f ( — (9 ■ curl Afc)(n • curl A°)—^(0 • curl A°)(n • curl Afc) j ds 

(3.25) - L f [ a }(9 • n)(iwAfc r + V T V fe ) • (AAf + V T lf) ds, 

icj j r 

where (B/~(0), £4(0)), (B^(0), £7 Z °(0)) are the material derivatives of (A k ,Vk) and 
(A^VJ 0 ) respectively. Now we will compute term by term (3.24). Remark at first 
that 


\ 0 (n d ) (a ?>), vf(0); Afc, -14) = 0 
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because the shape derivatives (A^(0), V®'(0)) vanish as the potentials (A^VJ 0 ) in 
the deposit-free configuration do not depend on Q d - This, together with ( |3.3| ), also 
implies 


B?(0) = (0 • V)A° + (VO) 4 A? and tf z °(0) = 6/ • W z °. 


Hence, by substituting (B^(0), ?7 Z O (0)) with the above expressions, one gets 

A(ffc) (Afc, V*; Bf(0), -iyM) = A(ffc) (Afc, U; (0 • V)A° + (V0)‘A?, -0 • Vy°) 

= <Si — ^2 


with 

S\ = f — curl Afc • curl ((0 • V)A° + (V0)‘A?) da; 
dOd A 

S 2 = ~ [ u(iwAfc + V14) • (iw((0 • V) + (V0) t )A? + V(0 • W, 0 )) dx. 
1UJ 

We compute S\ and S 2 


Si = f — curl Afc • curl(curl A° x 6) dx 

Jn d V 

= f curl(— curl Afc) • (curl A^ x 6) dx + f —(curl Afc x n) • (curl A^ x 6) ds 

Jn d M Jr V 

= f a(iujAk + VVfc) • (curl Aj 3 x 0) dx + f —(curl Afc x n) • (curl Aj 3 x 0) ds, 
Jn d Jr M 


and 


5 2 = 


— f <r(iwAfc + VVfc) • (iw(V(0 • A?) + curl A? x 0) + V(0 • VVf)) da; 
1W Jn d \ J 

= d f cr (iwAfc + VVfc) • ((iw curl A? x 0)+ V(0 • (iwA? + VVfc)) ) da; 

1W V J 

= — f cr(io;Afc + VVfc) • (icjcurlAj 3 x 0) dx. 

1CJ Jn d 


The last equality is obtained by integration by parts and by the fact that div (cr(ia;Afc + 
VVfc)) = 0 in Q d and that cr(icoAfc + VVfc) • n = 0 on T. Therefore 


A(ffc) (Afc, V k ; B?(0), - J7i°(0)) = Si - S 2 = l(curlA fe x ») • (curlA° x 0)ds 
(3.26) 

= J — ((0 • n)(curl Afc • curl A°) — (0 • curl Afc)(n • curl A°)^) ds. 
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Similarly, we have 

(3.27) 


Ao(a d )(A°,v) 0 ;B k (e),-u k (ej) 

= A 0 (n d )(Alv l °-,A r k (9),-vKe)) 


+Ao(n d ) (a?, Vj°; (6 • V)Afc + (V0)*A*, -0 ■ W fe ) 

= A o (a d )(A / fe (0),y fc '(0);Af,- : ^) 

+ J (^{9 • n)(curl Afc • curl A°) — (n • curl A k )(9 • curl A°) ) ds. 


From (3.24), (3.26) and (3.27), and considering the fact that the support of 6 is on a 
vicinity of T, we get (3.23). □ 

On T, we have 

curl A k • curl Aj 3 = (n • curl A k )(n • curl A^) + (curl A k x n) • (curl A^ x n). 

With the above equality and the relations ( |3.3| ), it follows that 

A Z' kl (n d )(6) 

= TJu [ |(1 - ^o)curlBfc • curlA° - 4(cr - cr°)(iwBfc + Vf/fc) • (iwA° + vv; 0 ) 
\ J \ Jn d IV V 

~YJy2 f {(7 “ 7o) curl (( 0 ' V)A fc + (V6>) 4 Afc) -curlA? 

\ J \ Jn d IV V 

-F(tr - <r°) (iw((0 • V)A fc + (V6») t Afc) + V(0 • Vt4)J • (kuA° + W°) J dx 


jw-n){ 


(n • curl A^)(n • curl A^) 


1CJ 

W 

"1 

V 

- — [a] (io; A kr + V r T4) • (i^A° r + V- 

1CJ 


( — curl Afc x n) • (^ curl A^ x n) 
V V 


v 0 )} 




ds. 


We follow the method of Hadamard representation to give an expression of z'ki(todW) 
dependent of ((A'(0), V r (Q)) or (B(0), U{6))) of the solution (A, F) by introducing 
the adjoint state (P i,Wi) G Q related to the solution (A^,VJ°) in the deposit-free 
case. The adjoint problem writes 


(3.28) 


5*(P Z , Wi; $) = L*(*, $) V(*, $) G Q, 


where for any (A, F), (\l/, 4>) in Q we have <S*(A, F; \F, <F) := <S(\F, <F; A, F) and 

L*(¥, $) := [ (1 - 4) curl A? • curl^ 

</o d V V 

+ [ hc7-a°)(ku A°+W°) • (icj\l/ + V4>) dx. 

•/fid 1CJ 

From the above considerations we easily derive the jumps condition for the adjoint 
states (P 1 ) 


[n • curlP/] = 0 

[/i -1 curlPj x n] = — ( — ^)curlA° 

V V 


on r, 


x n 


(3.29) 


on r. 
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It is worth noticing that the adjoint state Pj satisfies the Coulomb gauge condition. 
We are now in a position to express the results of proposition |3.5| with the use of 
the adjoints states (P i,Wi). Indeed, we have the following proposition (which is an 
immediate consequence of Proposition |3.5| and the definition of the adjoint state) 
Proposition 3.6. Let ( A 14) be the potentials induced by the coil k of the eddy- 
current problem with deposit domain Lid, (A®,V®) the potentials induced by the coil 
l for the deposit free case, and ( Pi,Wi ) the adjoint states related to which 

satify the adjoint problem (3.28). Then under the same assumptions as in Theorem 


3.2 for (i and a, the impedance shape derivative (3.28) can be expressed as 


A Z' kl {U d )[6) = ^ jjn-Q) 


(3.30) 


(n • curl Ak)(n • Pi — n • curl A®) 

— [/jl\ curlA/- x • ^^-(curl P^) + x n -^ curl A® x 

+ 3 [a](ic oA kT + V T V k ) • (iwP jT + WrWi + iujAl + V T Vj°) j da. 
1UJ J 


3.4. Explicit shape gradient formula. The shape derivative of f (Lid) is in 
the form 

(3.31) f'(^d)(0) = - 77 ^ [ (n ■ 9)gds, 

Kr Jr 0 

where the shape-dependent function g depends on the solutions to the forward problem 
(AfcjVfe), (A^, VJ°) and the adjoint state (P i,Wi). More precisely, g = gn + $21 for 
the absolute mode, and g = gn — #22 for the differential mode. For any l and k we 
have 


9kl = 


=f 




{ 


(n • curl Afe)(n • curlP/ — n • curl A^) 


— III I — curl A/c x n • —^ curlPj x n -- curl A, x n 

\T J \T° T° 

(3.32) + 2-[<x](iwA fcT + V T 14) • (iwP jT + V T W ; + iwA° T + V r y°) j) dC 


We choose the shape perturbation 6 such that 6 = gn on the interface T, which is a 
descent direction since 


f'(W) = -^| \d\ 2 ds < 0. 

4. Numerical algorithms for the deposit reconstruction. We recall that 
the computational domain Ll is a cylinder that contains the tube and the SP. We 
introduce a family of triangulation Th of Ll, the subscript h stands for the largest 
length of the edges in Th- The tetrahedrons of Th match on the interface between 
the conductive part (i.e. tube and SP a 7 ^ 0) and the insulator part (a = 0). The 
triangulation of the conductive parts (with deposits region) is given in Fig. [ 2 ] (see 
(A)-(B) for a real image and (A’)-(B’) for its F.E model). Since the variational space 
(of the regularized variational formulation) is based on H 1 functions, the numerical 
finite elements approximation will be based on nodal finite elements for the electric 
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Fig. 2. Pictures downloaded from @-Westinghouse: http://westinghousenuclear.com” . Figure 
(A) presents non clogging (healthy) SP, while Figure (B) presents a fully and partially clogged SP. 
Figures (A’)-(B’) are representations of used meshes for each configuration respectively. 


vector potential V as well as for the magnetic vector potential A. We shall mainly use 
Pi Lagrange nodal elements for both. In addition, the boundary conditions (A • n = 0 
on 30) are taken into account via penalization of degrees of freedoms that belong to 
30. Indeed the same numerical approximation procedure is applied to adjoint states. 

We now describe the gradient descent algorithm steps, the geometrical parametriza- 
tions and the procedure to accelerate iterations steps. The deposit is assumed to be 
located on the outer part of the tube and is concentrated (for the non axisymmet- 
ric examples) in one opening part of the SP (see Fig. [2]-(B)). The reconstruction 
is based on an intuitive approach, which consists in iteratively PO-approximating 
the geometry of the deposit on a predefined 3D grid. This method avoid to re¬ 
construct the mesh at each inversion iteration. The predefined grid is defined by 
•A4 = {T h d c Th, s.t V simplex K G with a facet parallel to Tube}, where h 
stands for the resolution of the grid. We give in Fig. [3]a clipping of A fh- We present in 



Fig. 3. Sketch of the deposit shape reconstruction (clipping of 3D representation) on a grid 
Gy. The reconstruction uses an invariant grid where PO-interpolation over this grid is set to cary 
out the new shape profile. 


































16 


3D shape reconstruction from eddy-current measurements 


Algorithm [l]the instances of an adapted step gradient descent. It is well known that 
the fixed step gradient descent algorithm converges if the step is sufficiently small. 
In our case we will allow the step descent to be large at least for the first iterations, 
and if the algorithm fails to maintain the decreasing of the cost functional, the step 
is reduced by a given factor. 1/2 < S < 1. The final geometry is the one for which no 
local variation (on the predefined grid) decreases the cost functional. 


Algorithm 1: Gradient descent algorithm 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 


Data: The impedance signal response of the tested configuration 
Result: Optimal shape approximation using interpolation on 3D grid 
Input: The resolution of the predefined grid h, ,Threshold: e 

n 1 

Input: T u Table of size N, a real t and — < 5 < 1 

2 

Input: P = E(Zog(0.4 * h)/log(8 )) 

Build the 3d mesh grid A : m x h X rriyh X N 
Evaluate the cost function f(T^) and the gradient Vf(L^) 
k = 0 

while || Vf (L k ) || 2 > e do 

t k = h/ ma Xl < n <tv |Vf(T^)| 
for 1 < p < P do 

tP = 8 X tP- 1 

LP = L k - tPVf(L fe ) 

Project the L,P (to the nearest value) on the predefined grid; Evaluate the cost function f (L,P) 
if f (L,P) < f(L k ) then 

Update f (L k + 1 ) = f(T p ) 

Evaluate the gradient Vf (L k ~^^) 

Break 

end 

if p=—P then 

I Print “A singular point is attained” 

I Exit () 
end 

end 

k = k -j* 1 

end 


5. Numerical implementation and validation. Numerical validation of the 
presented method is considered in this section. We use the software FreeFem++ m 
to deal with the finite elements discretization of the problem. We run our script on 
a cluster with distributed memory configuration. We use a direct matrix-inversion 
of the linear system where the factorization is achieved using sparse parallel solver 
(MUMPS HE]). We present and explain in the sequel some particular techniques to 
achieve performance of the direct eddy-current solver (and consequently the inverse 
solver). 

At each probe position we have to compute a solution associated to different 
source term. In order to (numerically) ensure divergence free condition for the source 
term one has to exactly mesh the support of the coil. If we build a new mesh related 
to the new probe position, we have to assemble new matrices and solve new systems, 
which are extremely memory-consuming. We therefore avoid this by creating and use 
a unique mesh that incorporates all possible probe positions in a scan of the tube. 
This allows us to only modify the right hand side of the system at each coil position. 
The factorization of the matrix is done only once per iteration. In order to further 
accelerate the resolution we also parallelize the matrix assembly since the cost of 
this part appeared to be the more expensive part if not done in parallel. Particular 
attention must be taken for the non-homogeneity (change of the conductivities and the 
permeability in the domain): We declare the variables a and /i as PO-Lagrange finite 
elements that depends on the elements labels of the non-partitioned mesh. Then, 
we apply a graph partitioning (e.g. scotch m or metis [19]) to create automatically 
partitioned new mesh. Since the partitioning process changes the elements labels to 
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|FA| : axisymmetric |F3| : axisymmetric 




Fig. 4. History of the impedance responses of the deposits during iterations in the axisymmetric 
configuration. \FA\ measurements (left) and |F3| measurements (right). 


the ranks of the used group of processors, we define the PO-Lagrange non-homogeneous 
domain variable on the non-partitioned mesh and then include them in the variational 
formulation that admits the partitioning (see m for more technical details). 

Numerical experiments deal with several configurations of test cases. Mainly we 
present an axisymmetric configuration, then we add the SP and consider the case 
where one of the SP foils (flow path) is clogged. 

The geometry of the computational domain includes a tube with respective inter¬ 
nal and external radius 9.84 mm and 11.11 mm. The coils are modeled by a crown 
with respective internal and external radius 7.83 mm and 8.50 mm. Both coils have 
length 2 mm and are separated by 0.5 mm. The scan step of coils is fixed to 1 mm 
and cover 20 positions along the tube, which length has been limited to 30 mm. 

We used the following values of the electromagnetic parameters. The frequency 
uo = 2007T, the magnetic permeability of the vacuum /iq = tt x 10 —6 , magnetic per¬ 
meability of the tube ft t = 1.01/io, the magnetic permeability of the SP ftsp = fto 
and the magnetic permeability of the deposit fid = The conductivity is taken 
<j t = 1 x 10 3 for the tube, ctsp = lx 10 2 for the SP and ag = cr t for the deposit. 

In all numerical experiments, the initialization of our algorithm takes a deposit 
with the lowest layers in the grid Gh i.e. with depth 0.5 mm equal to h : the precision 
of the fixed grid. 


5.1. Axisymmetric and non-axisymmetric geometries. In this part we 
consider two configurations of deposits in the vicinity of the tube: deposits around 
the tube far from SP and a deposit in one opening of water traffic lane of the SP. The 
first case, represents an axisymmetric configuration [18] and the second case represents 
a non-axisymmetric configuration because of the presence of SP and the deposit. We 
present in Figure [5] a slice on the plane (x,z) of the 3D computational domain. We 
show the shape of the axisymmetric deposits and the estimated deposits result 
of the inversion algorithm. Together with this plot we add the y-component of the 
solution Efc to show the penetration of the electromagnetic wave inside the tube and 
the deposits. With respect to &, a series of measured responses of the estimated 
deposit is presented in Figure [ 4 ] This shows the convergence of the method in the 
sense of minimizing the misfit function (2.3) presented in Figure [6] 

A more complex configuration consists in taking into account the presence of SP 
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Fig. 5. Slice on x-z-plan of the computational domain showing the plot of the y-component of 
the eddy-current solution and the deposits DJ (on the left) and the reconstruction in red (on the 
right). 



Fig. 6. Objective function with respect to iterations in the axisymmetric case. 


and therefore non symmetric deposit. The results for this configuration are presented 
as follows: In Figure [8] we plot a slice, on the plane (x,z), of the y-component of the 
solution E k together with the shape profile of the deposits and its estimation 0^ 
result of the inversion algorithm. The series of the impedance signal responses are 
given with respect to k in Figure [7] This highlights the convergence of our algorithm 
even with the presence of noise in the non symmetric solution (y-component of E&) 
as it can be seen in Figure [8] and also on the left plot of Figure [9j 


5.2. Arbitrary deposit shape. In this subsection in addition to the presence 
of the SP, we consider the reconstruction of a deposits with an arbitrary shape that 
does not match the parametrization used for the inverse problem: see Figure [lOj The 
results of the inversion algorithm is given (in terms of k ) in Figure[l3j The convergence 
in the sense of the impedance response measurements is given in Figure [TT] The 
minimization of the objective function with respect to the iterations of the inversion 
is presented in Figure [T2| 
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|FA| : with SP |F3| : with SP 




Fig. 7. History of the impedances during iterations for the SP configuration. \FA\ measure¬ 
ments (left) and |F3| measurements (right). 



Fig. 8. (left) Slice representation of the computational domain of the configuration with SP 
and a plot of the z-component of the eddy-current solution on x-y-plan. (right) Misfit function in 
terms of the inversion iterations. 


Appendix. Some useful differential identities. 

(A.la) curl (V/) = 0, 

(A.lb) div(curli;) = 0, 

(A.lc) (u • V)u = (V'u)'u, 

(A.Id) curlu x v = (Vu — (Vn) t )v, 

(A.le) V(u • v) = u x cuilv + v x curl u + (u • V)v + (v • V)u, 

(A.If) curl(u x v) = ndivv — vdivix + (v • V)u — (u ■ V)v. 


Appendix A. Proof of Lemma |3.1[ We develop the proof of the shape deriva¬ 
tive calculus presented at Lemma [XT] 
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Eel Y 

1 .0004562 
0.0004 
0.0003 
0.0002 
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6.2931e-8 



Eel Y 
0.0004562 
1 0.0004 
' 0.0003 
' 0.0002 
jO.OOOl 

6.2931e-8 



Fig. 9. Slice representation of the computational domain of the configuration with SP and a 
plot of the y-component of the eddy-current solution on x-z-plan. 




Fig. 10. Graph of the arbitrary shaped deposit that clogs one opening (foils) of the tube SP. 


Proof. By definition, one has 
A(n*) (A, V; V, $) *• / - curly A(S2 e ) ■ curly d y 

JQe V 

+ t - [ + V y V{tt e )) • (i^WN) dy. 

1U Jn„ 


With the variable substitution (Id + 6 ) 1 : y n i and the identities (3.1) related 
to A cur i, -4 c ]iv and VV ■ we rewrite the above form on a fixed reference domain Q = 
(Id + O)~ 1 Q ,0 as 


A(n e ) (A, V ;¥,$) = [ - (/ ) - curl A curi • curl fr cm , da; 

JoM |det(J + V0)| 

+ 7 - [ <71 det(/ + V0)|(I + ve)- 1 ^ + V0)“* (iwA curl + W v ) • (iw¥ cur i + V$ v ) d®. 


If (B (0),U(0)), (rj(0),x(0)) are respectively the material derivatives of (A, V) and 
then one can develop the above form with respect to 6 by considering the 
developments m 


(A.la) 
(A. lb) 


|det(/ + V0)| = l + div0 + o(6>), 
(/ + V6>) _1 =I-V0 + o(0). 
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|FA| : with SP noisy shape deposit 


|F3| : with SP noisy shape deposit 




Fig. 11. History of the impedances in the case or deposit with arbitrary shape: \FA\ measure¬ 
ment (left) and |F3| measurement (right). 



iteration of the inversion 

Fig. 12. Objective function with respect to iterations for the case of arbitrary shape reconstruc¬ 
tions. 


Since (A cur i(0), Uv(0)) = (A(U),U(U)), (^curi(O), ^>v(0) = the terms of order 

zero with respect to 6 in the development give exactly A(U)(A, V; \l/, <h), while the 
first order terms with respect to 0 yield 


A'(ft)(A, V; $) = A(fi)(B(0), U(0); $) + A(fi)(A, V; rj(0), x(0)) + h + Za, 

with = I — (—div0 + V0 + (V0)*) curl A • curl^dx, 

Jn l 1 

(A,2) 

J 2 = — f a{dwOI -VO- ( VO) t )(iujA + VV) ■ (kAlT+W) Ax. 

Jn 


We will rewrite the volume integrals Zi, Z 2 in terms of boundary integrals. Using the 
differential identities (A.l) and the fact that (A, V) satisfy the conditions (3.9), one 
verifies 


(-div 01 + VO + (VO)*) curl A = - curl((0 • V) A + (V6>)* A) + V(0 ■ curl A) + ^cr(iwA + VV) x 0. 
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Eel Y 
0.000447 
J0.0004 
•0.0003 
10.0002 
^ 0.0001 
1.0633e-7 


L 



(k=0) 




Eel Y 
0.000447 
[0.0004 


Eel Y 
0.000447 
[0.0004 


Eel Y 
0.000447 
[0.0004 



(k=12) 



(k=13) 


Fig. 13. Slice representation of the iterations of the shape reconstruction: case of an arbitrary 
shape clogging one opening of the quatrefoil SP. 


Hence 


!,=-/- curl((0 • V)A + (V6»)*A) • curl * da; + X n + X 22 , 

Jn P 

where In = / — V(0 • curl A) • curl T' dx and X 12 = / <x((iwA + W) x 0) • curl* 

Jn l 1 Jn 

By Stoke’s theorem, one has 

Xn = f — div ((0 • curl A) curl^) dx = f — (0 • curl A)(n • curl \l>) ds. 

Jn l 1 Jdn M 

By integration by parts (with use of differential identities A.l), we verify 

X 12 = — 7 — [ (jU&yvQI — V0)(iwA + X7V) + (0 • V)(iwA + W)} • (Iw*) da; 

1UJ Jn 

+ — cr(0 • n)(iu;A + W) • (io;^) ds. 

Jan 

Therefore 

Xi — — f — curl((0 • V)A + (V0)*A) • curl ^ dx 
Jn V 

- — [ cr{(div 61 — V0)(icoA + W) + (0 • V)(icjA + VH)} • (io;^) dx 
1UJ Jn 

(A.3) + / — (0 • curl A)(n • curl\l>) ds + — I a(0 • n)(iu;A + WV) • (icj\l/) ds. 

Jdn P 1LU Jdn 


dx. 
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Now we compute the term X 2 

I 2 = — [ cr(div 61 — V0 — (V6>) t )(iwA + W) • (W) dx + X 21 + X 22 , 

1Ld Jn 

with X 21 = A / adiv0(iwA + VV) • V$dx, 1 22 = — f ct(-V 0 - (V0)*)(iwA + VV) • V$ 
1( ^ Jn 1( ^ Jn 

By integration by parts, one obtains 

X 21 =J~ f cr(6 • n)(iwA + VV) • V$ds - — / iwa(0 • V)A • V$dx 

f <r£> 2 $(iwA + VV) -6dx- — f <t£> 2 V0 • V$dx. 
iw J n iu Jn 

Using integration by parts and the fact that div (<r(iu;A + VU)) = 0 obtained by 
applying the divergence operator to (3.9)i, one verifies 

x 22 = -^ f cr(V0) t (iwA + VV) • V<j>dx + d f X> 2 $(iwA + VV) • 6 dx. 

Jn ^ Jn 

From the differential identities ( |A.1| ), one deduces also that 

D 2 V6 + (ve)‘w = (6 • V) W + (V0) 4 VV = V(0 • VV). 

The above equalities yield 


dx. 


I2 —T21 + X 


^22 


(A.4) 


— [ a(div6I~ V0- (V0)‘)(iwA + VV) • (iw’J')dx + F [ cr(6 • n)(iwA + VV) • V$ds 

Jn Jan 

) 

- — [ criw((0- V)A + (V0)*A) • V$dx- 2. f crV(0 • VV) • Vl>dx. 

1( ^ Jn 1CJ Jn 


(A.3), ( |A.4| ) and the fact that cr(icoA + VU) • n = 0 on dQ imply 
h + x 2 =A(fl) ( - (0 • V)A - (V0) 4 A, -(0 • VV); 


(A.5) 


f — (0 • curl A)(n • curl \F) ds + f <j(Q • n)(icjA r + V r v) • (icj\l/ r + V<F r ) ds. 


From (A.2), (A.5) and the definition of shape derivatives (3.3), one concludes the 
result (3.10). □ 

Appendix B. Proof of Proposition |3.6[ We give the proof of the stated 
theorem 13.61 


Proof. Taking (\l>, <F) = (B/~(0), £4(0)) G Q in the adjoint problem (3.28) yields 
<S*(Pi, Wi; B fe (0), £4(0)) = L*(B k (0), U k {0)). 


On the other hand, taking (#,$) = (P i,Wi) in the variational formulation (3.18) for 
the material derivatives (B k (0),U k (0)) implies 


S(B k (0), U k (G ); P u Wi) = C(P h Wi). 
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Since 


<S*(Pi, Wi; B fc (0), £4(0)) = S(B fc (0), £4(0); P z , W) 

with the fact that div P/ = 0, one obtains 


L*(B k (e),U k (6))^C(Pi,Wi). 

= f — CUr l ((0 ■ V)Afe + (V0) t A/ c ) • curlP; dx 

Jn d V 

+ tj n <J ( iu ((^ v ) A ^(V») < A l )+V(e-W t )j -(iwPj + VWOdx 

(0 • n)(n • curl A/e) (n • curlPj) ds + A- f (0 • n) [cr] (icjA/~ r + V r T4) • (ic<;P/ r + V r Wi) ds. 

^ Jr 


L 


In f2\r one verifies 


(0 • V)A fe + (V0) t A fe = curl A& x 0 + V(0 • A/e), 
curl ((0 • V)Afc + (V0) t A/ c ) = curl (curl A^ x 0). 


Thus, considering (3.29)4 and (3.29)5, we compute 


L*(B fc (0),£4(0)) 

= Z+ [ — (0 • n)(n • curl A/e)(n • cnrlP/) ds# A- [ (0 • n)[a](iujA kT + V r T4) ■ (i^P/r + V r Wi) ds, 
Jr L^J 1cj Jr 

(B-l) 

with X = / — curl (curl A/e x 0) • curl P/dx + — / crkj( curl A/e x 0) • (icjP/ + VWj) dx. 


We remind that (curl A& x 0) belongs to X(fi). We multiply (3.29)i by (curl A/e x 0), 
integrate by parts and then take the complex conjugate, which implies 

X = j (— - ^r) curl A^ • curl(curl A/e x 0) dx — i f [cr](io;A^ +■ VVJ°) • (io;(curl A k x 0)) dx 

Jn d V V m Jn d 


L 


— curlP/ • 
r LM 



d s + [ 

"1" 

x n) 



Jr 



curl A^ • ((curl A/e x 0) x n) ds 


= f - curl A?- curl((0- V)A fc + (V6>)* A fe ) dx 
Jn d LmJ 

- - f [cr](iwA° + W,°) • (iw((0 • V)A fe + (V0) 4 A fc ) + V(0 • W fc ) ) dx 
1W Jn d V / 

(B.2) 

-J r (0-n){ 


— curl A/e x n ) • f -^r(curlP/)+ x n ) ds 

n J \n° I 























H. Haddar, Zixian Jiang and M.-K RIAHI 


25 


The last equality is due to the transmission conditions (3.29)2 - (3.29)3 for Pj and 
those for A& on T: [n • curl A] = [/j J ~ 1 n x curl A] =0. (B.l) and (B.2) imply 


■*(B k (0),U k (6))- [ - curl A?-curl ((0 • V)A fe + (V0) 4 A fc ) dx 

Jn d LM. 

- f [<r](iwA? + VVJ 0 ) • (iuUd ■ V)A fc + (V0)*A fc ) + V(0 • W fc ) ] da; 
1W Jn d \ ) 


+ — 


L (e ' n) {[k 


— (n • curl Afc)(n • curlPj) — [/j] ( — curl A& xn • —^ (curlPj) + x 


1 , 


n 


(B.3) 


+ — [<r](iu;A fer + V r T4) • (kjPj r + V r Wi) > ds. 


Considering the definition of L*(-, •), we substitute the above integral (B.3) in the 
expression of shape derivative of A Z^i (3.28) and finally obtain (3.30). □ 
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